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1 Introduction 

A series of elegant papers in the seventies has explored and clarified var- 
ious aspects of the large order behavior in perturbation theory in quantum 
mechanical and quantum field theoretic systems [|l], 0. In some cases this 
newly acquired knowledge was successfully used to obtain more accurate re- 
sults from the perturbative series, via, e.g., Borel summation method. In 
other cases, especially in four-dimensional field theoretic models, the ques- 
tion of what the sum of a badly-behaved perturbation series means, remains 
unanswered. For a review, see [Q. 

Recently, there has been a considerable interest in the so-called optimized 
linear delta expansion [4-17]. In simple quantum mechanical models, the 
method involves a rearrangement of the Hamiltonian as 

H = y + ^(?) (1-1) 

= Ho + 6H'\s=i, (1.2) 

where 

i^o-^ + yg^ H'^V{q)-^q\ (1.3) 

where f2 is a trial frequency. The model is then expanded in 5 up to a given 
order A^ according to the standard perturbation theory. The exact result for 
any physical quantity should not depend on the trial frequency Q artificially 
introduced above, but the corresponding finite order result S^ does. 

A proposal used often to fix Q is to use the "principle of minimum sen- 
sitivity" 0, i.e., to require that S^ be as little sensitive as possible to the 
variation of Q, 

If-"' t^'^) 

which determines Q and hence Sn order by order. Another possible criterion 
is the so-called "fastest apparent convergence" condition (see W^) 



For the anharmonic oscillator, V{q) = gq'^/4:, the optimized delta expansion 
with both prescriptions gives surprisingly good results already at the lowest 
nontrivial order, for any value of the coupling constant. 

The method seems very powerful and, in our opinion deserves a care- 
ful study. Already at first glance it displays several remarkable features. 
First, the method uses the standard perturbative technique. This makes the 
method potentially applicable to complicated systems such as field theoretic 



models of fundamental interactions (Quantum Chromodynamics, Quantum 
Electrodynamics, or the Weinberg-Salam theory), as well as to a wide class 
of quantum mechanical systems. 

Secondly, the method is nevertheless nonperturbative with respect to the 
usual coupling constant since the order-dependent determination of Q intro- 
duces a non-analytic dependence on it. The trial frequency Q, adjusted so 
as to eliminate the higher-order disaster of the standard perturbation the- 
ory, may be interpreted as a "vacuum parameter," similar to a certain field 
condensate or a physical mass parameter in the effective Lagrangian method. 



IS 



A particularly intriguing statement made by Duncan and Jones [^ 
that the optimized delta expansion (apparently) cures the problem of large 
order divergences even in non-Borel summable cases. If this were true in field 
theoretic models it could shed some important light on the instanton physics 
and vacuum structure of Quantum Chromodynamics. 

Finally, the method can be regarded as a generalization of the standard 
variational method, even though beyond the lowest order it involves no true 
variational principle. 

Despite these remarkable features, the theoretical basis of the success of 
optimized delta expansion remains unclear. A key observation may be that 
the method fails badly in a class of cases, where tunnelling effects are impor- 
tant (e.g., for the low- lying energy eigenvalues of the quantum mechanical 
double- well). If a method works in some cases but not in some other cases, 
it should be possible to understand the reason for that, and through such 
an analysis, to gain a better understanding of the mechanism underlying the 
success and of the limit of its validity. 



An important step forward has been made recently |T6[: It was demon- 
strated that the optimized delta expansion converges to the correct answer in 
the quantum mechanical anharmonic oscillator as well as in the double well, 
at finite temperatures. At zero temperature (hence for energy eigenvalues. 
Green functions, etc.), however, the convergence proof does not apply. The 



proof was then extended fl^ in a zero dimensional model (an ordinary inte- 
gral), to the logarithm of the integral which is an analogue of the connected 
generating functional. 

In this paper, we prove that the delta expansion for any energy eigenvalue 
of the quantum mechanical anharmonic oscillator. 



-" = TT + Try + TV ' (1-6) 



'1 ^ 

converges to the exact answer, if the order- dependent frequency Q is chosen 



to scale with the order N (asymptotically) as: 

xj^ = — = CN\ (1.7) 

where either 

^<7<^, C>Q, (1.8) 

or 

7 = 1/3, C > acg^^\ a^ ~ 0.5708751028937741. (1.9) 



In a zero-dimensional analogue model 
Zig,u;) = / — = exp 



- ( \^'q^ + \9q' 



;i.io) 



a similar result holds (see Appendix A|) but with the scaling index in a wider 
range 

J<7<^, (1.11) 

or 

7 = 1/4, C>a,g"\ (1.12) 

where in this case Uc is given by Eq. ( A. 25 ). 



Our convergence proof however does not apply to the double well case 

(CU2 < 0). 

Our work has some formal resemblance to the works by Buckley, Duncan, 
Jones and Bender |T3|, |IB|, |T^ and has in fact been inspired by them, but 
nevertheless differs considerably from theirs both in the method of analysis 
and in the results found (our result for the extreme case with a = ac however 
constitutes a generalization of their results). 

The principle of minimum sensitivity or the fastest apparent convergence 
criterion does not play any central role in the present work: the crucial idea 
of what we call scaled delta expansion is to scale appropriately the trial 
frequency with the order of expansion. In this respect our philosophy is close 
to a similar idea expressed by Bender, Duncan and Jones [0 , but we carry it 



to an extreme (relying solely on it), whereas their particular scaling behavior 
was suggested by the optimization procedure for Z{g,u). 

The ideas underlying our proof can be explained as follows. For dimen- 
sional reasons any energy eigenvalue of the anharmonic oscillator has the 
form, 

Eig,uj) = uE(^-^y (1.13) 



The standard (asymptotic) perturbation series then reads formally 

£;(^-*)=u;;^c„(4) , (1.14) 



n=0 -^ 



where the coefficients c„ have the known large order behavior at N -^ oo 
ii, 

^ ' ^ r{n + K + l/2). (1.15) 



The delta expansion for Eq. ( |1.6| ), i.e., a perturbative expansion in 

H = q +^q , (1.16) 



is equivalent to a substitution in Eq. (|1.14|) 0, [T^ 



^ ^ fi h + (5 ^^^^ j = ^x[l - 5/3(x)]l/^ 

^-'^•^, (1.17) 

where 

x = -, (3{x) = l-\, (1.18) 

followed by the expansion in S. A subsequent rearrangement of the series 
yields 



E^^'^ = to lim SNix) 



N^oc 






n=0 -^ 

N~n 



A parenthetical remark, useful for numerical analysis, is that the delta 
expansion for the double well case [u^ < in Eq. ( p..6| )] is given exactly by 
Eq. (|1.19| ) (in particular with the same coefficients c„ of the anharmonic 
oscillator), the only modification being 



u -^ \/ —up-, 

a; — >• a; = r2/V— cu^, 

/3(x) -^/9Dw(a;) = 1 + 1/x^ 



This assertion may at first sight appear puzzhng but is actually self-evident 
since the sign change of <J^ causes only trivial changes in the coefficients order 
by order of the expansion in H' [see Eq. ( [1.16| )]. 



Coming back to the anharmonic oscillator, S^[x) has the form of the stan- 
dard perturbation series with its coefficients modified by the f2-dependent 
factors A^^^[x). Except for Jv^ \x) they are positive definite numbers such 
that 

< A(f)(x) < 1 for n > 1, (1.20) 

(we assume a; > 1 throughout this paper). All of them [including Aq (x)] 
converges to unity as A^ — > oo if x is kept fixed, while they all vanish [except 
for Alq \x) which diverges] if x ^ oo with A^ fixed. Sn{x) can thus be 
regarded as a sort of regularized (or stabilized) perturbation series. The idea 
is then to play with the rates at which A^ and x are sent simultaneously to 
infinity so that the convergence is ensured. 

The convergence at the upper side of the sum (n = A^) requires that 

CNA'^i'\x) ~ ATl/x^^-i ^ 0, (1.21) 



this suggests x be scaled as x oc A^"^, 7 > 1/3. On the other hand, to 
guarantee the convergence of Aq (x) in the large A^ limit it is necessary that 

(3{xf = {l-l/xY -^0, (1-22) 



this implies that the scaling index 7 to be less than 1/2. The detailed study 
of the next section fully confirms these expectations. 



This paper is organized as follows. In Section ^T!] we prove rigorously that 
the sequence {Sn} converges to the exact answer as A^, x — > cxo according 
to the scaling in Eqs. ( [1.7| ) and ( |1.8| ). The convergence proof for the case 
7 = 1/3 is given in Section |2.2| . The uniformity of convergence, as well as the 
manner in which S^ diverges when the scaling index is outside the conver- 



gence domain are discussed in Section 2.3. These properties are compared to 



the numerical calculations of the delta expansion carried out to high orders 
(iV ~ 100) in Section 2A. In Section p.5| we reinterpret our result in terms 



of Seznec-Zinn- Justin's order dependent mappings. 

We conclude the paper by discussing several features of the scaled delta 
expansion and their possible generalizations in Section ^. Appendix A gives 



a convergence proof of the scaled delta expansion in the zero dimensional 



(ordinary integral) case; [Appendix B| discusses the properties of the roots of 



an equation used in Section pTl. 



2 Scaled delta expansion for the anharmonic oscillator 

In this section the scaled delta expansion for the anharmonic oscillator is 
discussed in detail. 



2.1 Proof of convergence for 1/3 < 7 < 1/2 



z-plane 



B 



1 i/p 



Figure 1: Cq contour in the complex 2;-plane. 



First we prove rigorously that the scaled delta expansion for the anhar- 
monic oscillator converges to the exact answer for K-th excited energy level 
{K = 0, 1, • • •), if the scaling of xn is chosen in the range Eq. (|1.8| ). Essential 
ingredients for the proof are: 

1) that an energy eigenvalue of anharmonic oscillator satisfies a once-subtracted 

, ImEig') 



dispersion relation |T^, 



E{g,uj) 



CqUJ + 



/■o , , Imi 

J~oo a' (a' 



TC J~Qo 



g'ig' - 9) 



6 



q f°° ,, ImE(X) 
= coUJ + ^ d\ y 2.1 

Tioj^ Jo A(A + g/iv^) 

where Cq = K + 1/2 and in the second hne dimensionless energy and 
couphng constant 

Ei\)=E{g',u)/u, \ = -g'/u\ (2.2) 

have been introduced. This relation follows from the analytic structure 
of the energy in the complex coupling plane, and from the known large 
A behavior E{X) ~ A^/^. By expanding Eq. ( p.l|) in powers of g/uj'^ 
one finds an expression for the perturbative coefficients c„ in terms of 
an integral involving \m.E{X) which was useful in the determination of 
their large order behavior ||^; 

2) the behavior of Imi?(A) at small negative coupling constant (small pos- 
itive A) determined by the tunnelling factor [|l|, 

Im E{\) = ^j^^^-""'"' exp (-- j [1 + 0(A)] , (2.3) 



for the i^-th energy level; and 

3) the positivity and boundedness of Im£'(A)/A^ [|T8| , 

00 >Im^(A)/A2 >0, for A > 0. (2.4) 

[This is a minor technical assumption. For the following proof, Im E[\)/\^ G 
L^[0, 00] is sufficient.] 

From now on uo will be set to unity, uo = 1. The A^-th order delta expan- 
sion approximant for E, Sn, can be readily constructed by the substitution 
rule Eq. (|1.17|) applied to Eq. (p.l|), followed by the Taylor expansion on S 
up to 6^ {6 = 1 at the end). This can be expressed compactly with the use 
of Cauchy's formula as: 

Sn{x) = Son{x) + Sin{x), 

So^ix) ^ Cox i ;^^^^^^(1 - PzY^', 

Jco Zm z — \ 



where 



^. X W(1-W7)V2 

(1 — wy^''^ + sw 



z-plane 




i/p 



Figure 2: Cqi contour in the complex z-plane. 



and 



s = 



9 



a;3/?A' 



P{x) ^ 1 



x^ 



(2.7) 



and the contour Cq is a small circle around the origin (Fig. 1). Because the 
integrand has no pole at z = 1, the integration contour may be deformed as 
in Fig. 2.Q Noting that the first term of the integrand (1 in the numerator) 
gives the exact energy in Eq. (|2.1|) , we get the following expression for the 
remainder 



Rn{x) = E{g,uj) - Sn{x) = Ron{x) + Rin{x), 
dz {l-Pzf/^ 



Ron{x) = CqX 



Coi Ini z 



N+li 



Rin{x) 



9 



'■°° ds 
I3x'^ Jo s 



^(^^1 



1)' 

dz F{(3z,s) 
27iiz^+Hz- 



(2.8) 
(2.9) 

:2.10) 



where the contour Cqi now encircles z = and z = 1 but no other singular- 
ities (see Fig. 2). Note that in Eq. ( p.lO| ) the integration variable has been 



^Note that this deformation does not encounter any obstruction due to the cut starting 
at z = 1/(3 > 1 or due to the poles of F{Pz, s): see below. 



changed from A to s in Eq. 



and 
ImE(A) 



ttA 



X=g/{x3f3s) 



(2.11) 



w-plane 



^-6 



% 

Cop X 



'•---^^ 



Figure 3: Deformation of Cop contour (w-plane) into cut contribution for 
Ron- 



Let us consider Ro]\r{x) and Rin^x) separately. First consider 



RoNix) = coxp 



N+l 



dw (1 - w)l/2 

Cop 27riw^+\w-p)'' 



(2.12) 



where the integration variable has been changed to w = (3z. The contour 
can now be deformed further to wrap around the cut on the real positive 
axis [1, oo) (see Fig. 3). The circle at infinity does not give any contribution. 
Introducing a real variable m by w = 1 + m ± ze, one finds 



Ron{x) 



CqxP 



N+l 



TT 



du 



u 



1/2 



Now 



-Roiv(a;)| < 



cqx(3 



N+l 



TC 



oo y 

(l+uW+^ 



(m + 1-/3)(1 + m)^+i' 



(2.13) 



/-1/2 



TT 



B{l/2,N+l/2), (2.14) 



9 



where 5(1/2, N + 1/2) is the Euler's Beta function. Using the Stirhng's 
formula we see: 

\Ron{x)\<^^^[1 + 0{1/N)]. (2.15) 



Next consider 

Rmix) 



9/3 



N 



ds 



x^ Jo s 



G(s)j^ 



dw F{w, s) 
Cop 27riw^+\w-p)^ 



(2.16) 



where w = j3z. In deforming further the integration contour the poles of the 
function F{w, s) [see Eq. (p^)] must be taken into account this time, as well 
as the cut. See Fig. 4. Accordingly Ri^^x) can be split as 



Rin{x) 



X JO s ,SZZ^„„ 



j=poles 

E i?S(a;)+i?S*H^ 

i=poles 



where 



7^, = —Res 



F[w,s) 



(2.17) 



(2.18) 



W=Wi 



is the residue of the i-th pole Wi, i.e., at the i-th root of the equation 

{l-wf^^ + sw = 0, (2.19) 

and 



C 



1 



du 



su 



1/2 



TT Jo (1 + m)^-i(1 -(3 + u)[u^ + s2(1 + uf] ■ 



^2.20) 



The bound on R]^^ '{x) can be set easily as follows. Obviously (1 — /3 

l/x^). 



\C\ < -lis), lis) ^ / 
n Jo 



du ■ 



su 



1/2 



But /(s) is continuous and behaves as 

71 

/(s) ^ — as s — i> 0, and /(s) 
o 



u-^ + s^il + u) 



TT 

— as s — i> oo, 
2s 



(2.21) 



(2.22) 



so that it is bounded in [0, oo). The contribution of the cut to -RiAr(x) is thus 
bounded by [see Eq. ( gTzD , Eq. (^1^)] 



R^in\x] 



< const. CigP 



N 



(2.23) 



10 










w-plane 



Ti-m 



1 



'•"•«.— i*^ 



o 



Figure 4: Deformation of Cq^ contour (ti;-plane) into cut and poles contri- 
butions for -RiAT. 



where we have used 



~ ds 



G{s) 



1 
n Jo 



°° ImE(A) 



Cl 



16 



{2K^ + 2K + 1), 



^2.24) 



(ci is the first coefficient of the standard perturbative series). 

To set the bound on J2i=poics ^ilfi^) requires a httle more work. The 



properties of the roots of Eq. ( p. 19 ) which he on the first Riemann sheet 



(they are the only ones that interest us) are studied in [Appendix B| . There 
are two such roots wi and W2 which move around on the first Riemann sheet 
as shown in Fig. 5 as s varies from to oo. At s = 0, they are coincident at 
Wi = W2 = 1; when < s < Sc = 3-\/3/2, they are complex conjugate to each 
other and lie on the two branches of the apple-like curve in Fig. 5; at s = Sc 
they coalesce again at Wi = W2 = —2. For larger values oi s, s > Sc, they are 
real and stay on the negative real axis: — cxd < wi < —2, and —2<W2< 0. 
Note that they keep always outside the integration contour Cop of Fig. 4. 

Near its i-th pole w ^ Wi the function F{w, s) behaves as 



F(w.s) 



2w^ 



1 



Wi + 2w — Wi 



(2.25) 



11 



[where Eq. ( |2.19 ) has been used] hence 



with 



7^i 



/H = 



Wi + 2' 



w^-\w-(3)' 



^2.26) 



^2.27) 



A httle obstacle arises in setting the bound on J2i=poies -^iM^^) due to the 



large residues near s = Sc where Wi ~ —2. Actually (see [Appendix B|) , in the 
vicinity of —2 the two roots are at symmetric positions with respect to —2, 
i.e., 

Wi,2 = -2t4v^/v^, 5 = 3^3/2 + 6, (2.28) 

so that these two large residues cancel in the sum. At this point it is con- 
venient to split the s integration range into three parts: (I) < s < Sq, (II) 
So < s < Sc, (III) Sc < s < oo, where sq = 2w2/3 is the value of s at which 

wi is purely imaginary wi{so) = \/3i (see Fig. 5). Any other choice of Sq < Sc 
would work. 




Figure 5: Behavior of the poles Wi and W2 in the complex w-plane as the 
parameter s varies. 
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(I) < s < So. 

In this region, two poles are complex conjugate to each other, W2 = w*, 
therefore 



E^' 



IN 



X 



i=l,2 



<2 



R^j^yx] 



(I) 



(I) 



For the right hand side, we find 



RlN\^, 



(I) 



< 



gp^ f ds 

X^ JO s 



fM 



Wi + 2 



ds 



Jo s 



< cig(3 



N 



^2.29) 



(2.30) 



where inequalities 



hi|>l, \wi-p\>l-p 



x^ 



\wi + 2\ >2, 



(2.31) 



(see Fig. 5) have been used. In deriving the last line of Eq. (|2.3CI| ), the 
integration region of s has been extended to [0, cxd). 

For regions (II) and (III), the difficulty of infinite residues at s = Sc can 
be overcome by rewriting the pole contribution to Eq. ( |2.17| ) as: 



ER 



(i) 
IN 



Xl 



i=l,2 



{II) + (III) 



2 4!r*GwlM^4^ + »^r^Gw/(-2)E ' 



i=l,2 ^ "'^0 S 

j: mix) + Rfiix) 

i=l,2 



Wi + 2 



X J sq S 



i^,2 «^i + 2 
(2.32) 



with obvious definitions for i?^*jv(x) and R[^{x). Bounds on each of Rilf{x) 
and R\]y{x) can now be set without difficulty. 

First consider i?j^(x). The function h{s) = J2i=i,2 Vl""^*!-^) + 2] is con- 
tinuous for s in [0, oo) (including s = Sc) and lim^^oo h{s) = 1/2, hence is 
bounded. Therefore 



i?S(x) 



gf3^ [^ ds ^, , cxgfi^ 

< const. _,, „ / — G(s) < const. — rr^r- (2.33) 

so s 2" x^ 



2%2 
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Next consider R^l^(x) {i = 1, 2) each of which is well-defined. We first 
note that 



l/H-/(-2)l 



dw' f'{w') 



<max|/'H|-Lr, (2.34) 



where F is a path connecting the points —2 and w and Lr is its length. Let 
us choose as Fj the trajectory w = Wi{s') traced by the i-th root as s' moves 
away from Sc = 3a/3/2 to the current value of s. Thus we define: 



Ti = <ww = Wi{s'), Sc < s' < s> , s > Sc, 

Ti = Iw W = Wi{s'), S < s' < Sc( , S < Sc 

(see Fig. 5). L^^ then satisfies 

Lr,(s) = \wi{s)+2\, s>Sc, 
Lviis) < A'\wi{s) +2\, s < Sc, 

where 



:2.35) 



A' = max 



LrAs) 



s£[so,Sc] \Wi{s) + 2| 

is some finite constant independent on s. It follows then that 



f{w,)-f{-2) 



Wi + 2 



< const, max !/'(«;) I 

^ i 



(2.36) 
(2.37) 

(2.38) 



in all cases. We use this relation and 



f\w) = -f{w) (^^LA + 1 



w 



w-(3 



(2.39) 



to get the bounds for w^l^{x) {i = 1, 2) in the regions (II) and (III). 
(II) So < s < Sc- 

In this region Wi and W2 are still complex conjugate, so that it suffices 
to bound R\p^{x) [see Eq. ( |2.29 )]. Using the fact that \wi\ is an increasing 
function of s (see [Appendix B| ) we get: 



thus from Eq. (p:39| ) 



kil > VS, \wi-j3\>Vl, 



,.„ 2(iV + l) 



(2.40) 

(2.41) 
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Use of Eq. (|23^ ) and Eq. (|0^ ) then yields 



R^j^yXj 



(11) 



< const. 



< const. 



3^/2x2 J«o S 



G{s) 



37V/2 



N 



X^ 



[2A2) 



(III) Sc < s < 00. 

In the region, a useful fact is that 

i/(»)i 



is a monotonically increasing function of w for negative w. For the smaller 
root, wi < —2, this implies that |/(wi)| < |/(— 2)|, hence 



m^ax|/'H|<|/(-2)|(^- + -]<-. 



^2.44) 



It follows that 



R^ikx) 



gNP^ foo ds 
(HI) < ^^^1 ~^^'^ 



< 






N 



For the second root W2 {—2 < W2 < 0) we instead use the fact that 



W2is) = --[l-Ms)?^'<--. 



Then 



yielding 



I/Ml < |/(-i/.)| 



max|/'H| < |/(-l/.) 
i 2 



2s 



N-l 



f3+l/s'' 



siN-l)+'- 



(2.45) 

(2.46) 
(2.47) 



< 



2s 



Af 



p y^^ps 



This leads to an inequality: 



^S(x) 



(III) 



< 



2^/3^-1 /-oof^s 



X' 






(3s 



^(--H' 



(2.48) 
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where 



In = 



J Sr S 



(2.49) 



and function G{s) is defined in Eq. ( p.ll| ). The contribution from large 
s is potentially dangerous. Fortunately, for s > Sc and at large x (as we 
are interested in the limit x -^ oo) A = g/{x^j3s) is always small, so the 
asymptotic estimate for Imi^(A) of Eq. ( |2.3|) can be used: 



/^2K+l 



K+?,/2 



x^(3\'^' 1-^ 



In 






- 



< 



v^7r3/2i^! V g 

^2K+1 ^3^ Af+i^+3/2 / 



Oig/{x'Ps)) 



x/27r3/2ir! V4 



N 



x^(3 



T{N + K + 3/2) [I + Oil/N)], (2.50) 



where K is the energy level. Substituting this into Eq. ( p. 481) and using 
Stirling's formula, we find 



Ri]^[Xj 



(HI) ttK\x^/3 \ 3gNj\Aex-^J ^^^^^ '^ 

(2.51) 

We are now in a position to set an upper bound to the full remainder Rj^ 
by using Eqs. (^, (|2T7| ), (^^ and ( |01 : 



Rn\ < -RoTV + 



p(cut) 



+ 2 



i?(^) 



lAf 



(I) 



+ 



5(0) 



+ 2 



5(1) 



(11) 



+ 



5(1) 



(III) 



+ 



5(2) 



(III) ' 

(2.52) 

with individual bounds given in Eqs. ( PJ[3|) , { ^^ , (j^]^), (^), (g^, 
d^) and (|3TD, 



|i?07v(x)|<^5SCill 






i?S(a:) 



R^j^yXj 



Ri^[Xj 



Ri^[Xj 



v/^iVV2 
< const. CigP^ 

< cig/3'', 

cigP"" 



Oil/N)] 



(I) 
< const 



(11) 



< const 



2^x2 ' 
cigN(3'' 



(2.53) 



(III) 



< 



cigN/3 

2^x 



3iV/2^2 

N 



2 ' 



USU^gN^ f Ax^\ f3gN\ 

(HI) 7ris:!a;2/3 I 3gN ) Uex^j 



TV 



[l + 0(l/iV)]. 
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If the trial frequency x = VL/uj is scaled as 

x = CN^ 
and the limit A^ ^ oo is taken at fixed K and g, then 



\Ra 



■OAf| 



(III) 



(2.54) 

3, (2.55) 

(2.56) 

if 7 > 1/3. This completes our proof for 1/3 < 7 < 1/2. 

Up to now no hypothesis was made on the constant C, which could be 



p(cut) 


7 


d(1) 


(I)' 


5(0) 


5 


5(1) 


(II)' 


5(1) 


while 












5(2) 


(III) 


^0 


) 





(2) 
IN 



X] 



(III) 



m 



chosen to depend on g and K. Actually the upper bound of \R 

Eq. ( |2.51| ) has a g^ dependence which may give a large remainder for large 

g. This suggests that we take 



a/3^7^ 



X = ag^'"i\ •, (2.57) 

to compensate such a dependence on g. We shall see in Section B^ this choice 



indeed gives a faster convergence. 

From the above results it follows that for the sequence {SN{xiy)} to con- 
verge to the exact energy eigenvalue it suffices that x^- = ^n/^ he anywhere 
within the range, 

C^p^l/3+e^ <xn< C2iV^/^~'^ (2.58) 

where the constants Ci,2, ^1,2 > are all fixed as A^ varies, so that there are 
actually an infinite set of sequences {Sn{xn)}, all of which converge to the 
exact answer. In the next subsection we shall refine the lower bound of the 
convergence range. 



2.2 Convergence for 7=1/3 

We analyze now the case with scaling index 7 = 1/3. Let us set 

X = ag^ 
To prove the convergence for 



U/3jYl/3_ 



it is sufficient to study 



Rl]y{Xj 



a > ar 



— G[s) 



(HI) x^ Js, s 



f{w2)-f{-2) 
W2 + 2 



(2.59) 
(2.60) 

(2.61) 
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since all other pieces of the remainder in Eq. ( p.53| ) are bounded by 

as N ^ oo. [R\/f{x)]^ui) can be estimated by the saddle point approximation. 



Note [as was done in deriving Eq. (|2.51|)] that at s > Sc and at x large, 



A = g/{x^(3s) is always small, so that the semi-classical estimate Eq. ( pl3|) 
can be apphed, yielding 

The saddle point in s is then determined by a compromise between the fac- 
tor \/w2 ^^ inside /(u'2) which tends to push it towards s = 00 (where 
W2 ~ ~'^/s —>■ 0), and the tunnelling factor exp[—Ax^l3s/{3g)] which tries 
to prevent it. With the scaling given in Eq. ( |2.59| ) the leading factor in the 
integrand for large A^ is 

c^{-lfexp[-N(^{a,W2)], (2.64) 

and 

0(«,^2) = ^^ ^^ + log|w;2|, (2.65) 

where we have set (3 = 1. The saddle point equation is then approximately 
given by 

^</'(«,^2)L„. = 0. (2.66) 

We studied this equation numerically, observing 1) that W2* is a decreasing 
function of a, and W2* -^ as a — > and W2* —>■ —2 as a — > 00; 2) that 
4>{a, it;2*) is an increasing function of a and 

(j){ac,W2,) = 0, (2.67) 

for 

ac = 0.5708751028937741 ■ • • . (2.68) 

(For a = ac, W2* = -0.2429640299735202 ■ ■ -, s* = 5.703557953364256 ■ ■ •.) 
The equations Eqs. ( p.66| ) and (|2.67| ) are essentially the same as the ones 



given m |^. 

Property 1) implies that the saddle point is always inside of the integra- 
tion region [recall that for Sc < s < 00, —2 < W2{s) < 0]. Property 2) shows 
that (f)(a, W2*) > for a > ac, hence 



^f^(x) 



^^^^^ ~ (3^ exp[-N<f){a, W2,)] ^0, if a> ac. (2.69) 



At a = ttc, convergence still holds due to the factor /5 in front of the integral, 



Ri^{x^ 



~ ,^JV ^ 



iVV3 



_ /3'^ ~ exp — — -. 0. (2.70) 

(III) ^ ^ I a2^2/3 ) V ; 



Furthermore {Sjy} converges if a = a{N) is a sequence such that liniAr^oo C({N) 
Ooo > and 

a{N) > ac, VN. (2.71) 

This result is of particular interest: in fact, the particular scaling index 
7 = 1/3 is known to follow from the fastest apparent convergence condition 
of the delta expansion for the energy [ffl. 



X 



a.g^/'^N^/^ + OiN-^/^), (2.72) 



as well as from the principle of minimal sensitivity for the partition function 
of anharmonic oscillator iTBH , 



x = a^ g^^^N^/^ + 0{N^). (2.73) 

Using the fact that the coefficients of the sub-leading terms in Eqs. ( |2.72|) 
and ( |2.73| ) are positive Q, |TB|, we thus complete the proof that the delta 
expansion converges to the exact energy eigenvalue for either of these choices 
for X. 

Summarizing the result of this subsection, the lower bound for the index 
for convergence is now weakened to xn > xnc where 



xnc = ac9 



i/3^i/3_ ^2.74) 



2.3 Uniformity of convergence of the scaled delta ex- 
pansion and divergence of R^ for scaling indices 
outside the convergence domain 

The uniform convergence in an arbitrary finite domain g ^ [0, G] follows if the 
upper bound for the remainder [Eq. ( p.52| ) and Eq. ( |2.53| )] is an increasing 
function of g. This indeed holds in general when C is an non- decreasing 
function of g, such as C = const, and C oc g^^^. 

Similarly, bounds on the remainder are increasing functions of K for g and 
A^ (> K) fixed. It follows that convergence is uniform for K < Kq < N. This 
behavior of the bounds suggests also that convergence is slower for higher 
energy levels. 
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To see that the scaled deha expansion actually diverges for 7 > 1/2, let 
us consider the first piece of the remainder, Rqn{x). We will find the upper 
bound in Eq. ( |2.15| ) is almost saturated for large N. Consider the difference 
of the upper bound and \Ron{x)\ [see Eq. (|2.14|) ]. 



CqxP 



N+l 



vr 



-B{1/2,N+1/2)~\Ron{x)\ 



co(3 



N+l 



nx 

N+l 



du ■ 



u 



-1/2 



< 



Co/? 



TTX 

N+l 



du 



M + 1 -/?)(! +m)^+i 

1 



mV2(m+1-/?) 



Co/3 



(2.75) 



The last quantity is negligible compared to the upper bound for 7 > 1/2 and 
for large A^. Therefore 

On the other hand, all other pieces in Eq. ( p. 53 ) are bounded by a constant 
[note when 7 > 1/2, [3{x)^ ^ 1 as A^ — > cx)]. This shows that the total 
remainder (hence Sjq itself) diverges as 

for 7 > 1/2 and N ^> 00. It is easy to see that the behavior Eq. (|2.77| ) also 
holds for 7 = 1/2. 

The divergence of the delta expansion for 7 < 1/3 on the other hand 
arises from \R{^{x)\(jii), since in this case the factor /?^ strongly suppresses 
all other terms of Rn in Eq. ( p.53| ). To show how it diverges we apply the 
semi-classical approximation Eq. ( p.63| ) and make a saddle point estimation 
of the s integral as was done in Section p.2| . 

The saddle point equation is now approximately given by 



d_ 
ds 



Ax'^(3s 
3^ 



A^logs 



0, 



^2.78) 



where the asymptotic form, W2 ~ — 1/s, s ^ 1, has been used. This is 
justified a posteriori from the resulting saddle point: 



t^ * 



Ax^j3' 



;2.79) 



which grows indefinitely for 7 < 1/3. Substituting this into Eq. ( p.61| ), we 
find that 

^e-^lJ^E)". (2.80) 

(III) V 4x3 y "- ' 



Ri^[x) 
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Thus the delta expansion diverges violently with alternative signs when N -^ 
cxD with 7 < 1/3. We recognize it as a heritage of the standard large order 
behavior of perturbative series (to which the delta expansion is reduced at 
X = 1), which is nicely tamed at a larger scaling index. 

Finally, at the critical index Eq. ( |2.59| ) but with a < ac, the discussion 
of the previous subsection shows that S^ diverges as 

c^i-lfexp[-Ncj){a,W2.)], (2.81) 

where 0(a,W2*) < for a < ac- 

2.4 Comparison with numerical calculation of ^^v 

To corroborate our proof and verify the general features found in previous 
sections, we have made a rather detailed numerical study of Sn- The method 
used is essentially the one described in [^] with some refinement. 

We first note that the A^-th order approximant Sn Eq- ( 1.19| ) can be 
written as 

N N-n / o \ n 

Sn{x)=xY.Y1 ^Kn 4 /^(^)'' (2.82) 



n=0 fc=0 ^^ '' 

where 

_ r(3n/2 + fc-l/2) 

''" ~ r(3n/2 - l/2)r(A; + 1) "■ ^ ^ ^ 

To compute the delta expansion for anharmonic oscillator [or for the dou- 
ble well potential, see the discussion below Eq. ( |1.19| )] to high orders therefore 



it suffices to find standard perturbative coefficients for the anharmonic oscil- 
lator Cn to the desired order A^ and use the relations above. 

We define the perturbative coefficients for the energy eigenvalue and the 
moments as: 

oo oo 

n=0 n=0 

The coefficients c„ and a^^^ then satisfy the following recursion formula [^ ^ : 



n + r 

^^-1) u-2) , 2(£+l )^ (,) 8(£ + 3) 
4(£ + 2) " ^ £ + 2 ^^0 

with the initial condition 



„('«) = iii,:i#a('-) + i^Ll±LY.Crat. - i^fif a,!S'.(2.85) 



Co = i\: + -, 4°^ = 1, a(f) = for n > 1. (2.86) 
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Eq. ( p. 85 ) has been derived by using the Feynman-Hellmann relation, 

(2.87) 



dE 
dg 




and the hypervirial theorem 



^i{f - l)(g^-2) + 2(£ + \)i,q\E - V)) - (g'+ V) = 0. (2.88) 

Note that this recursion formula is universal (i.e., independent of g): once 
c„'s hence ek,n^ are computed for each K^ they can be used repeatedly to 
study the delta expansion for various g and for different x. Our numerical 
results for the delta expansion approximants up to 100-th order have been 
produced in this way.0 Actually, with our approach, the coefficients tk.n up to 
500-th order were computed within 1 hour CPU time on VAXstation 4000/60 
(FORTRAN real*16 mode). 

Our illustrative results here will be limited to the case of the ground 
state energy of the anharmonic oscillator. For the numerical calculation of 
the exact energy, we used the method proposed in Ref. [0 and we computed 
the numerical value with an estimated relative error ~ 10^^^. 



s_gg (x) vs. X 



0.85 — 



0.80 — — 

5 10 15 20 

Figure 6: Typical plot of Sj^{x) [N = 99, (7 = 4). Exact energy (solid line) 
is also depicted for comparison. 



In |l6l a generalization of the recursion formula Eq. (2.85) to the delta expansion coef- 
ficients of E and (g^) is used, which, lacking the universality, is somewhat less convenient 
than our method from practical point of view. 
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A typical x-dependence of the delta expansion approximant Sn{x) at a 
fixed order is shown in Fig. 6 for N = 99, g = 4. An extremely flat region 
in a wide range of x, a very sharp rise at the lower end of the plateau, and 
a mild increase at large x, are clearly seen. These main features are well 
explained by the results of Section ^Tl] and Section p.3| . In particular the 
position of the sharp left end edge of the plateau is predicted by our analysis 
(Section ^) to be a^g^/^N^/^ which is ~ 4.19 for A^ = 99, ^ = 4. 
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Figure 7: a) the relative remainder cn = Rn/E with x = N^-^^ and g = 4: 
is plotted versus A^, b) same but log^g |eAf| is plotted. 



In Fig. 7a and Fig. 7b the relative remainder e^v = Rn/E with x = A^°' 
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is plotted versus A^ (in Fig. 7b log^Q |e| is plotted). They show that the 
scaled delta expansion (with an appropriate index) indeed converges rapidly 
to the exact answer. Fig. 8 illustrates moreover the general tendency that 
the convergence is faster, the smaller the scaling index is, as long as it is 
within the convergence domain, as suggested by the 7 dependence of the 
upper bounds in Eq. ( p.53| ). 

Fig. 9 shows a case with the scaling index below 1/3 {x = N^'"^). Sn is 
indeed seen to oscillate and diverge, as was expected. 

According to Eq. ( p. 771) , for scaling indices above 1/2, Sn should grow as 



^1 



X 



N 



Ari/2 



JV^-l/2^ 



(2.89) 



A logarithmic plot of eN (for three different values of 7, 7 = 1, 2, 3) versus 
A^ in Fig. 10 shows that such a simple law is indeed obeyed. 

The g dependence of Sn at a fixed order A^ = 40 is plotted in Fig. 11 for 
7 = 0.4. We examined two choices of C: C = 1 (dots) and C = {g/^Y^^ 
(solid line). There is a clear improvement in the second case as we argued 
at the end of Section ^]l] [Eq. (|2.57|) ]. It is also seen that convergence of 
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Figure 8: Plot of relative reminder ejv for fixed A^ = 40 and different values 
of the scaling index '-/ [C = 1, g = 4). 
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Figure 9: Relative reminder ejy versus N for x = N , g = 4: 
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ioqlO(abs(rei.err.)) vs N 
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— 



-2 ' ' ' ' ' — ' ' ' ' 

10 50 100 500 

Figure 10: Plot of log^^Q |e7v| (for three different values of 7, 7 = 1, 2, 3) 
versus N {C = 1, g = A). 



the scaled delta expansion is slowed at large g, as is suggested from the g 
dependence of the upper bound in Eq. (|2.53|) . 

All in all, our numerical study of the delta expansion fully confirms the 



results of the analyses in Section 01, Section |2]^ and Section p3|. 



2.5 Scaled delta expansion and order dependent map- 
pings 

The order dependent mapping |^ was proposed as a resummation method 
not based on the Borel transform. This approach involves a change of the 
expansion parameter g in Eq. ( |1.14| ) to a new variable A by a conformal 
transformation: 

g = pF{X), F(A)~ A + 0(A2). (2.90) 

The energy is then expanded in A instead of g: 



E = hm Sn{p), SMip) = /(A) Y: ^'Piip), (2.91) 



1=0 



with some prefactor function /(A); in the above Pi{p) is a l-th order poly- 
nomial of p with the order /. The parameter p in the mapping is then fixed 
order by order using some criterion (the exact energy should be p indepen- 
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Figure 11: Plot of log^g Mn versus g at fixed A^ = 40, 7 = 0.4. C = 1 (dots) 
and C = (^/4)i/3 (solid line). 



dent). In 0], the zero of Pn{p) with the largest module is chosen as p. This 
choice corresponds to the fastest apparent convergence criterion, Eq. ( |1.5| ). 

In the first paper on the order dependent mappings |Q, the equivalence 
of the linear delta expansion for the anharmonic oscillator (actually they 
called 9 instead of 6) and a particular order dependent mapping was already 
realized. In 0], they showed a mapping 



F(A) 



A 



;i-A)3/2^ 



with (in our notation) 



A = p{x) 



x^ 



P 



X[X^ 



[2m) 



(2.93) 



is equivalent to the delta expansion. To see this, note that the sum of the 
delta expansion up to A^, Sn{x) in Eq. ( |1.19| ) can be expressed in an another 
form (with uj = 1) 



N 



Sn(x) 



T.P'T.Cn 



T{n/2 + 1-1/2) 



VT^fo n=o 'T(3n/2-l/2)r(/-n + l) 



x{x'^ 



V. 



(2.94) 

A comparison of Sn{p) in Eq. ( |2.91| ) with Eq. ( p.94| ) shows the above equiv- 
alence. The only formal difference being the choice of the free parameter: p 
in the first case, and x in the second one. 
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Our proof of convergence of scaled delta expansion can thus easily be 
translated into the one for the order dependent mapping with a scaled pa- 
rameter p. It might be called a "scaled order dependent mapping." For the 
scaling x = CN"', one has p ~ g / {C^ N^"') for large A^. It follows from the 
proved range of convergence, Eq. (p..8|), that if p is scaled as p = gC /N'^ 
with 

1 < 7' < -, (2.95) 



or according to Eq. (|1.9D, as 



7' = 1, C'<\, (2.96) 



Q. 



then the order dependent mapping gives a sequence convergent to the exact 
answer (vice versa, a divergent one for 7' < 1 or 7' > 3/2). In particular, the 
extreme case with 7' = 1; C" = l/a^ corresponds to a particular choice made 
by Seznec and Zinn- Justin [^ (motivated by the fastest apparent convergence 
criterion) . 

For the zero dimensional case our convergence proof ( [Appendix A| ) implies 
that an order dependent mapping based on 

F{\) = jY^2^ (2-97) 

converges to the exact answer if p is scaled as p = gC /N"' with 

1 < 7' < 2, (2.98) 

or 

7' = 1 and C < \. (2.99) 

This follows from a relation A = I3{x) and p = g/[x^{x'^ ~ 1)] which holds in 
this case (see H)- 



3 Discussion 

Apart from those features specific to the delta expansion, the convergence 
proof presented in Section |2T2 relies essentially only on dispersion relation 
for the energy, Eq. ( p.l|) , and on the small coupling constant behavior of the 
imaginary part of energy Eq. ( |2.3| ) (determined from a semi classical estimate 
of tunnelling amplitudes). These are the same ingredients used in the study 
of the large order behavior itself pi. We therefore believe that it should be 
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possible to generalize our proof to a wide class of quantum mechanical as well 
as field-theoretic models. The way in which the trial frequency in introduced, 
however, might be different case by case. 

The lower bound for the scaling index 1/3 in the case of quantum me- 
chanical anharmonic oscillator reflects the physical dimension of the coupling 
constant, hence the dimensionality of the system D itself. A simple minded 
generalization of the proof would yield the lower bound 1/(4 — D) for the 
index, suggesting that the cases of higher dimensional theories {D > 2) re- 
quire a nontrivial extension of the method (e.g., modifying the way Q is 
introduced). 

As compared to the conventional approach using the Borel resummation, 
the present method is more direct, requiring no analytic continuation of the 
Borel transform. However it is still to be proved that all (or some subset of) 
cases in which the perturbative series is Borel summable can be treated by 
scaled delta expansion. 

As regards the interesting cases with degenerate classical minima where 
the standard perturbation series is not Borel summable, our preliminary 
numerical as well as analytical study on the quantum mechanical double well 
has not yet yielded a definitive answer on the applicability of scaled delta 
expansion. The convergence proof of Sections |2]l| and |2.2| clearly does not 
apply to that case as it stands. The main difference being 

/3dw = 1 + ^ > 1 (3.1) 

now. We hope to come back to this problem in a near future. 

Although the scaled delta expansion is quite different in spirit from the 
"optimized" delta expansion in which the trial frequency Q is determined 
order by order, for instance, by principle of minimum sensitivity, the results 
of this paper in no way diminish the virtue of the latter approach. Rather, we 
believe that our results put the optimized delta expansion on a much firmer 
ground. 

The remarkable (empirical) success of principle of minimum sensitivity 
for the anharmonic oscillator energy eigenvalues at lowest orders, might well 
be related to the proven existence of an infinite set of sequences that satisfy 

a, g^/^N^/^ <xn< C2N^'^-'^ , (3.2) 

all of them converging to the correct answer. This may be of practical impor- 
tance because in more complicated systems optimization procedure might be 
essential in getting a rapidly converging answer. 
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Also, as a by-product of our analysis, optimized delta expansion based on 
the fastest apparent convergence criterion (with xjsf lying just on the lower 
bound of Eq. ( |3l^ ) 0), has been given a rigorous proof of convergence. 

Let us compare also the results found here with those of Bender, Duncan 
and Jones []T^. In the zero dimensional case, they find that if x is scaled 



as N^/'^ — corresponding to principle of minimum sensitivity — then the 
delta expansion converges for both signs of uP. (See also |^.) We proved in 
Appendix A| that Sn converges with 7 in a much wider range (1/4 < 7 < 1/2) 



with an arbitrary proportional constant or 7 = 1/4 and C > acQ^^^. On the 
other hand it is not obvious in our approach to see the convergence for tu^ < 
with 7 = 1/4. Evidently, the two approaches are somewhat complementary 
in this regard. 

In the quantum mechanical case, the convergence index 1/3 (with a par- 
ticular proportionally constant) found by Duncan and Jones ||TB[ for the full 
generating functional at finite temperatures, is again just on the boundary 
of the convergence domain Eq. (|3.2D for energy eigenvalues. Therefore with 



their particular scaling [see Eq. ( |2.73| )], delta expansion is proved now to 
converge to the correct answer for energy eigenvalues as well. 

To conclude, there is still much to be clarified but the scaled delta ex- 
pansion appears quite promising as a new resummation method for the per- 
turbative series. Generalization of our proof to the anharmonic oscillator 
Green's functions, as well as to some other simple cases, will be discussed in 
a forthcoming paper. 
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Appendix A Convergence in the zero dimensional case 

In this appendix we give a convergence proof for the scaled delta expan- 
sion applied to the integral. 
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The delta expansion is defined by 
dq 



where 



exp 



2 ^ 



M-(.^ 



^')e 



-^gq 



N 



5=1 



N ^(fc)^o) 
fc=0 



{fmn^T.'-iT^5'- 



k\ 



(A.2) 



(A.3) 



As in the one dimensional case, the delta expansion is equivalent to the 
substitution 

u -^ y^l]2 + 5(cu2-fi2)^ g^5-g, (A.4) 

into the exact expression Eq. ( |A.1|) , followed by an expansion in 5. By using 
the dispersion relation for Z, 



Z{g,uj) = -Z{-g/uj 

UJ 



z(,,,) = irf^'-^w 



cj Jo TT A + g/uj^ 



(A.5) 



(where use was made of the large g behavior Z(g,uj) ~ g ^/^), the same 
procedure as in the one dimensional case then yields 



>Af 



where s 



pN+i f^ds ~ r dw 1 

7TUJX Jo s Tco 2iTi w^~^^ {w — (3) 

g/{ijj^x^[3X) and 



Fo(w,s), 



(A.6) 



Fq{w,s) 



(1 



w 



1 3/2 



[1 — wY + sw' 



X 



Vt ^ 1 



u 



x^ 



(A.7) 



and Cq is a small circle around the origin. This is quite analogous to S'itv 
considered in Section p.l| , and the proof is consequently very similar to that 
case. We deform the contour so as to wrap around the pole at w = /3, poles 
of the function Fo{w, s) and the cut running from w = 1 io w = oo. The 
residue of the pole at (3 gives exactly (the minus of) Z{g,ui) itself. (We shall 
set u; = 1 from now on.) The remainder i?jv = 2' — Zjy is given by the sum 
of the contribution of the poles of the function Fq{w, s) and that of the cut, 



R 



N 



R 



(cut) 

N 



p(cut) 



+ Rn 

f^N+1 



(poles) 



X 



where 

°° ds\m.Z r°° du 

s 



TC 



1 



1 



U 



3/2 



R 



(poles) 

N 



TT {1+U)^+^U+1/X^U^ + S{1+U 

1 1(1- w^y/^ 



pN+1 ^oo dslmZ ^ 

X Jo S T^ i wf Wi — (3 Wi + 1 



From Eq. ( A. 8 ), we get 



p(cut) 



<P 



7V+1 



(A.9) 



(A.IO) 
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where only the factors u + l/x^ and v?' have been retained in the denominator 
of Eq. (|A.8|) and use was made of the relation /g°° ds Im Z /{ns) = Cq = 1, Cq 
being the zeroth order coefficient of the perturbative expansion. 

The poles of the function Fq{w, s) are at: 



wi,2 = -(2-s=F v^ 



4s 



(A.ll) 



When < s < Sc = 4, the poles Wi and W2 are complex conjugate to each 
other and move on the unit circle centered at the origin, \wi\ = 1. When 
■5c < s < oo, they are on the negative real axis and — oo < Wi < —1, 
— 1 < W2 < 0. To study R^° , it is convenient to divide the integration 
over s in three parts. 

(I) < s < So, where sq (< Sc) is arbitrary. In what follows, a choice 
So = 2 is made for definiteness. In this range. 



E 



1 



1 (i-WiY/' 



-^ wf Wi- P Wi + 1 



< 2^l^x\ 



(A. 12) 



where use was made of \wi\ = 1, Itfj — /9| > 1 — /3 = 1/x^, |1 — Wj| < v^ and 
|w7i + 1| > \pl (note that Wi = ±i for s = 2). The contribution of this part 
to Eq. (|A.9|) is thus bounded by 



p(poles) 



i-N 



(I) 



< 2'/^x(3''+\ 



Next to avoid the large residues near s 
integrand of Eq. ( |A.9| ) as: 



(w 



(A. 13) 
— 1) we split the 



E 



f{wi 



Wi + l 



< E max If 



E 



Wi + 1 



-L 7, 



^/2A 

Wi + l\ ' 1 + /3 



+ 



(A. 14) 



where f{w) 



w 



ll/2, 



/[w {w — (3)], Lr- is the arc length between the 



points —1 and Wt, and we have used the fact that | J2i ^/{wi + 1)| < A {A 
being a finite constant). Thus the contribution of the last term in Eq. ( A.14| ) 
to Eq. ( [A.9| ) is bounded by 



R^''> 



i-N 



< 



V2A(3 



N+l 



x{l + (3) 



The contribution from the first two terms of Eq. ( |A.14| ), Rj^ 
bounded as follows. 



(poles) 



(A.15) 



can be 
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(II) So < s < Sc- We note LrJ\wi + 1\ < A' with a finite constant A' and 
1 |l-w;2|^/2 / jY I 1/2 \ 



l/'l < 



\W2 



N 



\W2-P\ \\W2\ \W2-P\ \W2 



< ^(^ + ^ + 271) 



(A.16) 



since \w2\ = 1, V2 < |1 — 102] < 2 and \w2 — P\ > I. Therefore 

/ 1 \ f}N+l 

< 2V2A' {n + 1 + ^] ^ . 

(II) V 2a/2/ X 



H 



(poles) 

N 



(A.17) 



(III) Sc < s < 00. Note that Lp. = \wi + 1| in this region. Also, |/| is a 
monotonically increasing function for Wi < 0. It follows that: 



l/'l < I/I 



A^ 



+ 



wi\ \wi-f3\ 



+ 



1/2 



< 



Wi\ 



V2 
1 + 13 



N 



1 V 

l + /?^4, 



(A.18) 

where we have used inequalities |/| < |/(— 1)| = -\/2/(l + /?), \wi\ > 1, 
\wi — P\ > I + P and |1 — Wi| > 2. Therefore, the contribution of Wi to 
Eq. ( |A.9|) is bounded as: 



R 



(III) 



<^'-D^ 



On the other hand, we have for W2 



\W2 



N 



+ 



1/2 



W2-l3\ \\W2\ |tf2 - /9| 1^2-1 



< 



V2s^ 



(A.19) 



."'^H. 



(A.20) 
where use was made of \w2\ > 1/s (following from W2 = —(1 — W2)'^/s < 
— 1/s), and 1 < |1 — W2\ < 2, \w — P\ > j3. Since in (III) the argument of 
Im Z is always small we may use the approximation. 



ImZ(A)~-^exp(-^)[l + 0(A)], A « 1. 

This leads to a bound for the contribution of W2: 
1 fAg^^ 



R 



< - , 

(III) vrx \x 



^ivr(iv + i)+ 1 + i r(iv) 



(A.21) 



(A.22) 



From the collection of the upper bounds, Eqs. (pTTSD , ( [AA5| ), (|AT7| ), (|A:T9| ) 
and (|A.22|) , we conclude that 



R 



'N 



(A.23) 
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if 

A^ ^ oo, x = CN^ -^ oo, with 1/4 < 7 < 1/2. (A.24) 

(2) 
Finally, for 7 = 1/4, a saddle point estimate of [i?)^ ]iii analogous to 

Section |2.2| shows that Rn ^ if 

a>ac = 0.972632810758477 • • • , (A.25) 

where x = ag^^^N^^^. On the other hand, the principle of minimum sensitiv- 
ity criterion [^ or the fastest apparent convergence criterion 0, |15[ requires 



a = 1.072985504616992 ■ ■ ■ > a^- (A.26) 

For this choice of a therefore we get an exponential decrease of the remainder, 
in accord with what has been shown earlier 11, P3. 



Appendix B Behavior of the poles of F{w,s] 



In this Appendix the behavior of the roots of Eq. ( p. 19 ) 



{l-wf/^ + sw = (B.l) 

which lie on the first Riemann sheet will be studied. A graphical analysis of 

{l-w)^ = sW (B.2) 

shows easily that Eq. ( [B.2| ) has either one (for < s < Sc = 3-\/3/2) or three 
(for Sc < s < 00) real roots. One real root which is always present and lies 
in [0, 1] (which we call w^ below) however does not satisfy Eq. ( |B.1| ). More 
precisely, W3 is a solution of Eq. ( p.l|) lying on the second sheet of the cut 
w-plane. Eq. ( p.l| ) has then either two real solutions (for Sc < s < 00) or 
two complex poles (for < s < Sc = 3v^/2) in the first sheet, which will be 
denoted as Wi, W2. 

At s = 0, Wi = W2 = 1 obviously. Near that point the behavior of Wi^2 is 

^^2 = 1 - s2/3eT2W3_ (B_3) 

The behavior of wi 2 near the confluent point, —2, is important in the 
analysis of Section 01. By setting 



s = 3^3/2 + e, wi,2 = -2 + 5, (B.4) 

in Eq. (|B.2|) and expanding up to the second order of 6 it is easy to find 



S = T^v^. (B.5) 
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The two roots thus stay at the symmetric points with respect to —2 in the 
vicinity of —2. The corresponding residues of the function F{w, s) thus 
behaves as 

9 V^ 

ResF ~±^, wi,2 2, (B.6) 

hence their sum vanishes when s -^ Sn- 



The behavior at s — ^ oo can be easily read off from Eq. (KT): a negative 



solution of Eq. ( |B.1|) should go either to — oo (wi) or to (^2)- Their leading 
behaviors aX s -^ 00 are: 

wi{s)-^~s^[l + 0{l/s)], W2{s)^--[l + 0{l/s)]. (B.7) 

From Eq. ( |B.2| ) we see that wi^2 satisfy 

W1W2WZ = 1, < W3 < 1, (B.8) 

from which we get an interesting property 

W1W2 > 1. (B.9) 

When wi^2 are complex conjugate to each other (0 < s < Sc) this means 

\wi\ = \W2\ > 1, (B.IO) 

which is also obvious from Fig. 5. 



Actually we can go a little further. From its implicit definition Eq. ( B.2 ) 
W3 is seen to satisfy 

dw-i 2swn ,^ 

OS 3(1 — wsj^ + 2s^W3 

It then follows from Eq. (P3.8D that the product W1W2 (or |wip in the complex 
region) is an increasing function of s, a fact also used in Section O. 
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